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Abstract Pulsating flows with a total reverse flow phase are ubiquitous in phys- 
iological systems in normal and pathological conditions. Irregularity of hemody- 
namic parameters in such flows is correlated with the appearance and develop- 
ment of several arterial pathologies. We study the relations between flow wave- 
form parameters and the wall shear stress (WSS) related quantities such as mean, 
root-mean-square, gradient of WSS and the oscillating shear index. The phase- 
averaged velocity profiles measured by the digital particle image velocimetry are 
used to estimate WSS utilizing the Womersley pulsating flow model. In addition 
to the Reynolds and Womersley numbers, another dimensionless parameter, pul- 
sating index (PI) which is the ratio of forward flow rate to the reverse flow rate is 
required. PI is essential for the complete description of the flow patterns with the 
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total flow reversal. We demonstrate significant effects on the WSS quantities due 
to the pulsating frequency and PI. Furthermore, the particle residence time (PRT), 
measured using the three-dimensional particle tracking velocimetry, shows signifi- 
cant correlation with WSS related quantities. For example, shorter residence time 
correlates strongly with the high pulsatility index. Correlations of other quantities 
are qualitatively similar though weaker, suggesting cross-coupling effects. It is pro- 
posed that PRT, estimated in-vivo from the trajectories of magnetic or contrast 
tracers, can provide an alternative diagnostics of wall shear stress parameters. 

Keywords reverse phase ■ retrograde fiow ■ pulsatile fiow ■ wall shear stress • 
particle residence time ■ digital particle image velocimetry ■ three-dimensional 
particle tracking velocimetry 

1 Introduction 

Flow patterns with a total reverse phase, i.e. a phase with a negative flow rate 
typically referred to as retrograde or regurgitate flows, are present in the cardio- 
vascular system. During flow reversal, the fluid flow direction is altered and the 
velocity gradients change sign abruptly. Local flow reversal occurs when the fluid 
in the near wall region is reversed, while the total flow rate of the cross section is 
positive. Aorta flow experiences a total reverse phase that results from the rapid 
deceleration of blood at the end of systole before the closing of the aortic valve [6l 
I13II15] . Increase in the amplitude or duration of the reverse flow phase is associ- 
ated with rapid blood flow elevation which can result from invasive procedure or 
arterial pathology [41112]. Haddad et al. [7] predicted a different form of local flow 
reversal, the so called "off wall reversal" , where the negative velocity is found at a 
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certain distance from the wall, while the near wall velocity is positive. One of the 
main interests in pulsating flows with the total reverse phase is the shear stresses 
developed on the arterial wall. Arterial wall shear stresses (WSS) are known as a 
main regulator of endothelial cell function and vascular structure [11114] . Irregular 
values of WSS properties are strongly related to cardiovascular pathologies. De- 
spite the clinical importance, a detailed analysis of the important hemodynamics 
parameters, along with the reliable measurements of WSS in flows with the to- 
tal reverse phase, is lacking. Although a variety of non intrusive imaging based 
measurements such as MRI, Ultrasound and particle velocimetry are available, 
accurate measurement very close to the wall are non-trivial to obtain due to the 
limited spatial resolution, light reflections and deficit of tracer particles in the 
area. Due to the high uncertainty in the near wall region data, the curve-fitting 
or model-fitting approaches are commonly used [lOj . 

Pulsating cardiovascular fiows with local fiow reverse are well described using 
two dimensionless parameters: the Reynolds and the Womersley numbers. For 
example, Finol et al. [3 found that the maximal values of non-dimensional mean 
wall shear stress and non-dimensional wall shear stress gradient increase with 
the Reynolds number, while Ref. [8_ demonstrated the significant effect on WSS 
magnitude of the Womersley number. We show that in case of the pulsating flows 
with the total reverse phase, a complete description of the fiow is not possible 
without an additional dimensionless parameter, the pulsatility index (PI), [5]. 
This index is an arterial blood-flow velocity parameter quantifying the pulsatility 
or oscillations of the fiow waveform, e.g. (T7] . 

We estimate the wall shear stress r, its root-mean-square (rmag), time gradient 
(WSSG), and oscillating shear index (OSI), using digital particle image velocime- 
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try (DPIV), and reveal their dimensionless relations. In addition, we measure 
Lagrangian trajectories of flow tracers using the three-dimensional particle track- 
ing velocimetry (3D-PTV) and demonstrate insightful correlations between the 
directly estimated particle residence time (PRT) and WSS for different flow con- 
ditions. In regions of recirculation (regions of increased plague deposits), particles 
have significantly longer times of residence and therefore this parameter was cor- 
related with hemodynamic properties in flows with a total reverse phase. 

2 Materials and Methods 

2.1 Experimental method 

A custom-design flow test rig was developed to create well-controlled pulsatile 
flows, as shown in Fig. [1] 

[Fig. 1 about here.] 

The pulsating flow is created in an elastic tube (made of Tygon B-44-4X, Saint 
Gobain, inner diameter of 3/4" and wall thickness of 1/8", L/d > 40, elasticity 
modulus of 12 MPa, estimated distensibility 1.36 x lO"'' Pa^^) by a set of three 
DC voltage driven computer controlled gear pumps. In order to allow high quality 
optical measurements, the tube was placed in a 800 x 300 x 200 mm glass tank 
flUed with a refractive index matched liquid (60%u, glycerin- water solution). A 
4:0%w glycerin-water solution was used as blood mimicking fluid. The flow was 
monitored using pressure transducers (EW-68075-02, Cole Parmer) and magnetic 
flow meter (MAG 1100, Danfoss). We created a set of flow waveforms (dipolar and 
unipolar sine w/o total flow reversals, unipolar sine, cardiovascular-like waveform, 
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see Fig.[2|) in the following range of the flow parameters: Reynolds number between 
60 900, Womersley number a : 6 -^ 11 and pulsatility index (PI) 1.5 9. 

[Fig. 2 about here.] 

We implemented DPIV using a dual Nd:YAG laser (532 nm, 120 mj/pulse, 
Solo 120XT, New Wave Research), a high-resolution CCD camera (12 bit, 4008 
x2672 pixels, TSI Inc.) as shown in Fig. [3^. Silver-coated hollow glass spheres 
(14 ^im, 1.05 g/cm^, TSI Inc.,) were used as seeding particles. More than 100 
images were acquired at each of the 7 phases of the pulse period, triggered by the 
pressure transducer. Phase averaged velocity was estimated using a commercial 
software (Insight 9.1, TSI Inc.) and verified with an open source software [16] 
(www.openpiv.net) . 

Lagrangian trajectories of fluorescent microspheres (nominal diameter 60 /im, 
p = 1.51 g/cm'^) were measured using 3D-PTV, as presented in Fig.[3)D, illuminated 
by the LED hght source {MB - LL306 - i? - 24, MetaPhase Tech.) and recorded 
using 4 high-speed cameras (8 bit, 1280 x 1024 pixels, Mikrotron GmbH) and 
stored on the digital video recording system (10 Industries) at rates of 100 -h- 250 
fps. Open source software (www.openptv.net)is used for image processing [11| . 
Detailed description of the experimental method is given in [9l[6] . 

[Fig. 3 about here.] 

[Table 1 about here.] 

2.2 Wall shear stress related properties 

Wall shear stresses were estimated using two methods: a) by derivation of measured 
near wall velocity using 3 points curve fitting, and b) by utilizing the Womersley 
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solution for pulsating flow in a rigid tube [19_. The flow rate waveform Q has been 
calculated by integrating the measured velocity profiles. The Fourier coefficients 
of the signal were calculated using fast Fourier transform (fc = 32, TV = 100): 



AT-l 



Qk= ^"^ 



n=0 

where Qn is the n*'* component of the flow rate waveform in its Fourier represen- 
tation. Velocity profiles are estimated from the Fourier representation of the flow 
rate waveform according to the inverse Womersley relations: 

uir,t)=Rely i ^ ^ i e^--' k (2) 

[ -^"^ \ J^/^a„Jo (j^/^Q„j - 2,/i (j3/2Q„j ; J 

where Jo Ji are the zero and first order Bessel functions respectively, uj is the 
pulsation frequency, j = \/--T and an is the n*'' Womersley number defined as 
an = R\/^ and R is the inner radius of the pipe. Finally, WSS can be estimated 
deriving the obtained Womersley velocity profiles at different phases. 

The fiow waveforms were characterized using three dimensional numbers: the 
maximal velocity Reynolds number -Remax: 

i?e = i?emax = C/max-D/;^, (3) 

where Umax is the maximal time averaged velocity, D is the tube inner diameter 
and u is the kinematic viscosity of the fiuid; the Womersley number a: 

a = R^JuJv (4) 

and the pulsation index (PI): 

PI = (Qmax — Qmin) /Qmean (5) 

Flow waveforms with total reverse flow phases are characterized by high PI val- 
ues. The following WSS related properties are considered in the present study: 
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dimensionless mean wall shear stress, 




Hereinafter the superscript * denotes the dimensionless parameters, normalized 
using the time period of the pulsation, T, and the equivalent Poiseuille wall shear 
stress Ts'. 



rs = (7) 

where U is the mean value of the time averaged velocity profile in the cross section, 
U = Qmean/TT-R^. Another hemodynamic parameters is the wall shear stresses 
magnitude, averaged over a cycle: 

r5iag = ^ ^^y|r(i?,i)|dtj , (8) 

where | • | denotes the absolute values of the WSS t{R, t). Next are the oscillatory 
shear index (OSI): 

OSI=ifl-^) (9) 
^ V ^mag / 

and the time averaged wall shear stresses gradient (WSSG): 

T 



WSSG 



-sT\tJ 



dT{R,t) 



dt 



dt (10) 



In addition to the measurements, we analyze the same properties for a larger 
variation of the dimensionless parameters Re, a, PI, using a simulation based on 
the Womersley solution, namely Re= (70, 170, 200, 350, 690); a= (6, 8, 9, 11, 13); 
PI= (1.5, 2.5, 4, 5, 6.5, 9). In each simulation run, only one parameter was changed 
while the other two were kept constant. Experimentally obtained WSS and related 
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quantities are compared in the following to the trends of these hemodynamic 
parameters estimated from the simulation. 

3D-PTV experimental method provides the velocity and accelerations of La- 
grangian tracers in the region of interest and during the time of pulsation. More- 
over, we deduce from the Lagrangian dataset the particle residence time (PRT), 
defined here as the time interval required for a particle to exit a spherical domain 
of pre-defined radius v o 

PRT{xo,t,ro) =t:\x{t)-x{to)>ro\ ■.x = {r,t) (11) 

We estimate PRT for a number of values in the range of t-q = 0-^0.25 mm, averaging 
PRT values of all measured trajectories in a given experimental run, and present 
in the following a comparative study of the flow patterns. 

3 Results 

3.1 Velocity profiles 

DPIV measured velocity profiles (symbols) and their corresponding Womersley 
solutions (curves), according to Eq.Q, are presented in figure [3] for 16 different 
runs and different phases, 6. Profiles at the beginning of the cycle and at the 
acceleration phases are rectangular in shape, as described in the literature for mid 
frequency laminar fiows [191118] . Near wall local flow reversal can be seen in most 
runs for 9 > 180°. Total flow reversal appeared in runs 11, 14, 16. It can be seen that 
on overall, the measured velocity profiles fit well with the Womersley solution. On 
average, the deviation of the measured velocity values from the analytical solution 
is below 18%. As it can be expected, the strongest deviations are found in the 



Wall shear stress in total reverse flow 



9 



near-wall region. WSS estimated from measured velocity profiles shows a good 
compatibility to WSS calculated from the inverse Womersley solution within the 
average uncertainty limits of 20%. 

[Fig. 4 about here.] 



3.2 Dimensional analysis of WSS properties 

Comparison between simulated WSS values and quantities calculated from mea- 
surements can be seen in figure O Dimensionless WSS related quantities are not 
influenced by the variation of Reynolds number if the Womersley number and pul- 
satility index are kept constant (not shown here). Mean WSS (rmean), according 
to the Womersley model, is constant and equal to 0.75rs (dashed line in Fig. [5^). 
Normalized WSS magnitude, shown in Fig. [5]d, maintains a constant value for PI 
= 1.5. For higher PI values, WSS magnitude increases with an increase in PI or 
frequency, when the effect of increasing PI is more significant. Oscillating shear 
index (OSI) is shown in Fig.[5j:. Minimal values are received for the lowest pulsatil- 
ity index tested, i.e. PI = 1.5 and are frequency dependent. OSI values increases 
substantially for the range of 1.5 < PI < 6, showing dependency of both PI and 
the Womersley number. For PI > 6, the OSI values increased rate is reduced and 
apparently saturates to an asymptotic value. WSSG maintains a quasi-constant 
value and it does not dependent on PI for the lowest range of Womersley numbers 
tested, (Fig. [5}1). However, for the higher frequencies, WSSG values increase al- 
most linearly with both frequency and PI. Measured WSS magnitude and WSSG 
usually follows the simulation trend, except for runs with substantial reverse flow 
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phase, where the discrepancies can not be explained by the measurement uncer- 
tainty. OSI values shows grater discrepancies compared to other quantities. 

[Fig. 5 about here.] 

3.3 Particle residence time 

Mean particle residence time, estimated according to Ea. (|lip . as a function of the 
sphere radius ro is presented in figure [6l It can be seen that particle residence time 
increases with ro, as expected, though at different rates. Comparing the residence 
time of two runs at the same Reynolds number, we observe shorter residence time 
for higher PI. The same trend is observed when we compare the runs at the same 
frequency - the residence time is shorter for higher PI. According to the analysis 
above, WSS related quantities are strongly influenced by pulsatility index (PI). 
This fact indicates a possible relation between the particle residence time and WSS 
related properties. 

[Fig. 6 about here.] 

In order to test this hypothesis, we plot the PI, OSI and WSSG versus PRT 
for all the runs in Fig. [7^-c at constant rp (200 /im). The results show very strong 
correlation between particle residence time and pulsatility index. Thus, runs with 
shorter residence time are characterized by higher PI values. The correlation be- 
tween the particle residence time and OSI or WSSG is not as good, nevertheless, 
the general trend is clearly seen - the particles residence time is short when the 
OSI and WSSG are high. 

[Fig. 7 about here.] 
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4 Summary and discussion 

DPIV measured velocity profiles show good agreement with the Womersley's an- 
alytical solution for pulsating flow in rigid tubes, with slight deviations of the 
measurements from the analytical solution near the wall. We can address the de- 
viations to the relatively lower quality of measured data due to light reflections 
and inherent inaccuracy of defining the wall location in DPIV images. Thus com- 
paring measured data to model profiles and applying image processing methods, 
the closest distance from the wall, for which PIV velocity measurement can be 
considered as reliable, is estimated as 0.7 mm or 0.07y/R. Therefore, the Womers- 
ley rigid tube flow model can be considered a suitable reference for measured flow 
properties. 

Using computational simulations based on the Womersley solution and the 
experimental results, we show that properly normalized WSS properties do not 
depend on the Reynolds number at the same frequency of pulsations. However, 
different values of WSS related properties were found for flows with an equal 
Reynolds number and Womersley number but for different waveforms. We con- 
clude that the two parameters are not sufficient to describe all the relevant flow 
properties and an additional parameter is introduced, namely the pulsatility index 
(PI) , calculated from the properties of the flow rate waveform [10] . Experimental 
runs with total flow reversal phase are characterized by high PI values. Simulations 
based on the Womersley model show that the WSS related quantities (except the 
mean normalized rmeari)i highly influenced by the two parameters: PI and the 
Womersley number. WSS magnitude and WSSG both increase with increase of 
Pl/frequency. PI is stronger correlated with the WSS magnitude, while the Wom- 
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ersley number affects stronger the WSSG. WSS magnitude and WSSG follow the 
trends obtained from the simulations, except for the experimental runs with total 
flow reversal. Measured OSI values are significantly lower than predicted values. 
OSI values may be influenced by the duration of the reverse flow phase - parameter 
that was not considered in the present study. 

Since direct WSS measurements in in-vitro and in-vivo applications are diffi- 
cult, the indirect non-intrusive methods, such as particle residence time (PRT) 
can become valuable measures of the local WSS related properties. PRT can be 
measured by the 3D-PTV method in flows with local or total flow reversal as 
the time required for a particle to exit the near wall region. Shorter PRT was 
found for the flow patterns with the highest PI, as compared to the runs with 
higher Reynolds/Womersley numbers and lower PI ratios. We revealed strong anti- 
correlation between PI and PRT, i.e. if a given flow pattern is characterized by a 
high pulsatility index, the measured particle residence time is expected to be low. 
Correlations between PRT and other WSS related properties are not as obvious, 
yet the trends of decreasing values with increasing PRT are clearly seen. Further 
study is needed, since it seems likely that quantitative correlation between PRT 
and WSS related properties could be established given a larger data base and a 
higher quality measurements in the near wall region. 

To summarize, the dimensional, empirical and computational analysis of WSS 
related quantities demonstrate that the use of a single Reynolds number together 
with the Womersley number is not sufficient to characterize the pulsating flows 
with total flow reverse. In the present study, a pulsatility index that takes into 
account the amplitudes of the forward to reverse flow phases was utilized. We 
realize that in some case of the pulsatile flows, the definition of PI should be 
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extended in order to include the time intervals of forward/reverse flow phases. In 
a view of a possible in-vivo estimate of the particle residence time, which is shown 
here to be strongly correlated with the WSS properties of interest, an additional 
research in this direction can be important. 
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Fig. 1 Hydraulics system schematics; the flow is created by the controlled gear pumps and the 
reverse flow phase is due to the additional gear pump installed in the opposite flow direction. 
Voltage control of the pumps and the feedback from the flow meter and pressure transducers 
are regulated using the PC 
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Fig. 2 Different types of generated flow patterns: (A) Unipolar sine with a secondary flow 
rate peak. (B) Dipolar sine. (C) Unipolar sine with total flow reversal. 
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(a) (b) 
Fig. 3 Velocity measurement systems: (a) PIV setup scheme (b) 3D-PTV experimental setup. 
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Fig. 4 DPIV measured velocity profiles (symbols) and the corresponding Womersley solution 
(lines) for 16 runs and different phases (color online). 
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Fig. 5 Normalized WSS related quantities variation with frequency and pulsatility index, (a) 
Mean WSS. (b) WSS magnitude, (c) Oscillatory index, (d) WSS gradient. Values calculated 
from the simulation are presented by lines (a =6 - blue solid line, a = 9 - red dashed line, 
a = 11 - green dotted line and o = 13 - black dash-dot-dash line) . Values calculated from 
direct measurements are presented with symbols, while triangular markers represent the sets 
with reverse flow phase. 
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Fig. 6 Mean particle residence time (PRT) for the three representative runs as a function of 
the interrogation sphere radius, rg. 
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Fig. 7 Particle residence time as a function of the following parameters: (a) pulsatility index 
(PI) (b) oscillatory shear index (OSI) (c) wall shear stress gradient (WSSG). 
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Run 


Waveform 


T (sec) 


Kemax 


a 


Pi 


Method 


3 


Dipolar sine 


2 


788 


9 


2.4 


PIV 


4 


Unipolar sine 


1 


806 


13 


1.6 


PIV 


5 


Unipolar sine 


4 


607 


6 


3.3 


PIV 


6 


Unipolar sine 


2 


540 


9 


2.6 


PIV 


7 


Square wave 


1 


630 


13 


1.7 


PIV 


8 


unipolar sine 


2 


942 


9 


1.5 


PIV 


9 


unipolar sine 


4 


173 


6 


4 


PIV, 3D-PTV 


10 


unipolar sine 


4 


213 


6 


3.9 


PIV, 3D-PTV 


11 


unipolar sine 


4 


60 


6 


8.9 


PIV, 3D-PTV 


12 


unipolar sine 


2 


353 


8 


2.4 


PIV 


13 


unipolar sine 


2 


192 


8 


3.2 


PIV,3D-PTV 


14 


unipolar sine 


2 


69 


8 


6.4 


PIV, 3D-PTV 


15 


unipolar sine 


1.2 


336 


11 


1.8 


PIV 


16 


unipolar sine 


1.2 


64 


11 


4.8 


PIV,3D-PTV 



Table 1 Summary of experimental runs for the various wave forms (unipolar, bipolar sine, 
square wave), periods, Remax,oi and PI. Runs 1 and 2 were performed at steady flow rate for 
the calibration purposes and are not shown here. 



